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Accelerator neutrino oscillation experiments seek to make precision measurements of the neu- 
trino flavor oscillations v ^ — ^ v e in order to determine the mass hierarchy of neutrinos and to 
search for CP violation in neutrino oscillations. These experiments are currently performed with 
beams of muon neutrinos at energies near 1 GeV where the charged-current quasi-elastic interac- 
tions vin — > £~p and i>ep — > £ + n dominate the signal reactions. We examine the difference between 
the quasi-elastic cross-sections for muon and electron neutrinos and anti-neutrinos and estimate the 
uncertainties on these differences. 



I. INTRODUCTION 

Since the invention of neutrino beams at accelera- 
tors and the consequent discovery of the two flavors of 
neutrinos [H, the reactions vgn — > £~p and Dip —> £ + n, 
which are the dominant reactions of muon and electron 
neutrinos with energies from 200 MeV to 2 GeV, have 
played an important role in studies of neutrino flavor. 
These charged-current quasi-elastic (CCQE) interactions 
are important not only because they identify the flavor 
of the neutrino through the production of the lepton in 
the final state, but also because the two body kinematics 
permit a measurement of the neutrino energy from only 
the observation of the final state lepton. 

Accelerator neutrino experiments like T2K0, 

NOvAQ and a number of proposed experiments seek to 

make precision measurements of the neutrino flavor os- 

•u *■ ( " ) ( " ) (-)(-). , . , . 
cillations v ^— > v 6 or v e — > v ^ m order to determine 

the mass hierarchy of neutrinos and to search for CP vi- 
olation in neutrino oscillations. Uncertainties on differ- 
ences between these cross-sections for muon and electron 
neutrinos will contribute to experimental uncertainties in 
these flavor oscillation measurements. 

Interactions of the charged-current with fundamental 
fcrmions, such as vgd — > £~u, have no uncertainties in 
the differences between the reactions for muon and elec- 
tron neutrino interactions because the weak interaction 
is experimentally observed to be flavor universal. In par- 
ticular, the effect of the final state lepton mass in this 
two body reaction of fundamental fermions can be un- 
ambiguously calculated. 

One such calculable difference occurs because of radia- 
tive corrections to the tree-level CCQE process. Radia- 
tive corrections from a particle of mass m in a process 
with momentum transfer Q are of order — log — , which 
implies a significant difference due to the mass of the final 
state leptonQ. Although this effect is calculable, it is not 
accounted for in neutrino interaction generators used in 
recent analysis of experimental data, such as GENIE 
NEUT @, If and NUANCE ^ . 



There are, however, cross-section differences due to 
lepton mass which cannot be calculated from first prin- 
ciples with current theoretical tools. The presence of the 
target quarks inside a strongly bound nucleon lead to a 
series of a priori unknown form factors in the nucleon 
level description of the reaction, e.g., vpn — > £~p. It is 
the uncertainties on these form factors combined with 
the alteration of the kinematics due to lepton mass that 
leads to uncertainties, and that is the focus of the results 
of this paper. 

There is also a modification of the reaction cross- 
sections due to the effects of the nucleus in which the 
target nucleons are bound. The model incorporated in 
GENIE [6], NEUT0J1 and NUANCE ^ is a relativis- 
tic Fermi gas modeljlOl Ull ] which provides a distribution 
of nucleon kinematics inside the nucleus. A more so- 
phisticated description from a nuclear spectral function 
model[I2] is implemented in the NuWro generator [l3j. 
We do not consider the effect of the nucleus in this work, 
although it may be important in the relative weighting 
of nucleon kinematics at low energy. However, this work 
provides a blueprint for studying the effect of the final 
state lepton mass in different nuclear models. 

II. NUCLEON FORM FACTORS 

The cross section for quasi-elastic scattering of neutri- 
nos at energies relevant for oscillation experiments may 
be calculated from the Fermi theory of weak interactions 
with the effective Lagrangian, 

£ c ff = —j= yJ^ x J^ H ^ + Hermitian conjugatej , (1) 

where Gf is the Fermi constant and the J are the leptonic 
and hadronic currents. The form of the leptonic current 
is specified by the theory to be 

J(t)\ = ^7a(1-75)V' I /« ) ( 2 ) 

because the leptons are fundamental fermions. However, 
as mentioned above the hadronic current for quasi-elastic 
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scattering depends on unknown form factors of the nu- 
cleons. The hadronic current can be decomposed into 
vector and axial components, 



tA _ tA 

J (H) — J V 



(3) 



Jy contains three terms related to the vector form factors 
Fy, F v and F v , and J a contains three terms related to 
the axial form factors Fa, F\ and F„. A description the 
the bilinear covariant structure of the currents is given 
in several standard texts and review papers [T34T6j . We 
follow most closely the notation of Ref. YWL 

From the effective Lagrangian of Eq. [T] and currents 
above in Eqs. [2] and O the quasi-elastic cross section on 



free nucleons is: 
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where the invariant Q 2 — — q 2 and q is the four momen- 
tum transfer from the leptonic to hadronic system, M is 
the mass of the nucleon, 9 C is the Cabibbo angle, and E v 
is the neutrino energy in the lab. The combination of 
Mandelstam invariants s and u can be written as, 



s- u = AME V -Q 2 -m 2 



(5) 



where m is the mass of the final state lepton. The func- 
tions A(Q 2 ), B(Q 2 ) and C(Q 2 ) depend on the nucleon 
form factors and £, the difference between the anoma- 
lous magnetic moment of the proton and the neutron: 
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Note that the form factors themselves are functions of 
Q 2 in Eqs. ME 

F v and F v are the vector and Fa and Fp the axial 
form factors of the first class currents. First class currents 
conserve both time and charge symmetry. In addition, 
first class vector currents commute with the G-parity op- 
erator while first class axial currents anti-commute with 
it [lij] • The terms associated with F v and Fa are consid- 
ered the leading terms in the hadron current since they 
have no dependence on the four-momentum transfer, ex- 
cepting that of the form factors, and they are not sup- 
pressed by powers of the final state lepton mass as Fp 
is. 

Vector elastic form factors are precisely known at 
Q 2 = from the nucleon electric charges and magnetic 
moments and are precisely measured over a wide range of 
Q 2 in charged-lepton elastic scattering from protons and 
deuterium. The axial nucleon form factor at zero Q 2 is 
precisely measured in studies of neutron beta decay. At 
higher Q 2 , much of our knowledge of the axial form fac- 
tors comes from muon neutrino quasi-elastic scattering 
measurements. For Q 2 ~ 1 (GeV/c) 2 , the vector form 
factors and the axial form factors are observed to follow 



a dipole form, that is 



F{Q 2 ) oc F(0)/(1 + Q 2 /C 2 f 



(9) 



where the constant C is often expressed as a mass of the 
same order of magnitude as the mass of the nucleon. At 
high Q 2 the vector form factors do not follow the dipole 
structure [18]. The neutrino scattering data contains con- 
flicting results among different experiments [l9l - [23| which 
limit our ability to effectively use that information to 
constrain the axial form factor. Pion electroproduction 
experiments [U Hfl have also measured the axial form 
factor at Q 2 < 0.2 (GeV/c) 2 . 

The form factor Fp is determined from PCAC which, 
under minimal assumptions, states that [26]: 



5^Ja = C4> 



(10) 



where <f> is th e renormalized field operator that creates 
the 7r + and C is a constant which may be computed at 
Q 2 = 0. PCAC gives the following relation between Fp 
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and the pion nucleon form factor, g-^NN, 



F P (Q 2 ) 



2M 2 F A (0) 



( F A (Q 2 ) 9*nn{Q 2 



\ F A (0) g„NN(0) (l + ^ r ) / 



,(11) 



where Mj is the pion mass. The Goldberger-Treiman 
relation[23] predicts 



9*nn(Q 2 )F„ = F A (Q 2 )M, 



(12) 



where is the pion decay constant. Under the assump- 
tion that the Goldberger-Treiman relation holds for all 
values of Q 2 , then Fp is 



Fp(Q 2 



2M 2 F A (Q 2 ) 
M 2 + Q 2 ■ 



(13) 



This is the relationship that is used in all modern neu- 
trino generators [(^M liH] . 

Fy and F\ are form factors associated with the second 
class current (SCC). The existence of such currents re- 
quires charge or time symmetry violation, and current 
measurements show the size of these violations to be 
small. Additionally a nonzero F v term would violate 
conservation of the vector current (CVC). Both F v (0) 
and F^(0) can be limited experimentally in studies of 
beta decay. Almost all current analyses of neutrino data 
assume that the SCCs are zero. The vector SCCs only en- 
ter the cross-section in terms suppressed by m 2 /M 2 , but 
there are unsuppressed terms involving the axial SCC 
form factor. 



III. MUON AND ELECTRON NEUTRINO 
QUASI-ELASTIC CROSS SECTION 
DIFFERENCES 

In this section, we will study the dependence of the 
muon-electron cross-section differences as a function of 
E v and Q 2 . Differences arise due to the variation of kine- 
matic limits due to the final state lepton mass, different 
radiative corrections to the tree level process and uncer- 
tainties in nucleon form factors. Equations |6] and [7] con- 
tain explicitly the dependence of the CCQE cross-section 
in terms of the form factors. Lepton mass, m, enters in 
both A(Q 2 ) and B(Q 2 ) and these terms involve all the 
form factors discussed above. Note that Fp and F v only 
appear in terms multiplied by m 2 /M 2 and therefore are 
negligible in the electron neutrino cross-section, but not 
in the muon neutrino cross-section. 

As metrics, we define the fractional differences between 
the muon and electron neutrino CCQE cross-sections: 



S{Eu ,Q 2 ) = f— 
A(£„) = 



JdQ 2 ^. 



JdQ 2 £p 



(14) 
(15) 
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FIG. 1. The total charged-current quasi-elastic cross-section 
difference for neutrinos (top) and anti-neutrinos (bottom) due 
to the kinematic limits in Q 2 . This difference is —A defined 
in Eq. 1151 meaning that the electron neutrino cross-section is 
larger than the muon neutrino cross-section. 

The integrals over Q 2 in Eqs. [T4l and [T51 are taken within 
the kinematic limits of each process, and those limits 
depend on lepton mass as discussed in the next section. 

Another useful metric is the difference between a cross- 
section in a model with a varied assumption from that 
of a reference model. Our reference model derives F v 
and F v from the electric and magnetic vector Sachs form 
factors which follow the dipole form of Eq. [9] with C — 
c 2 M v = (0.84) (GeV/c) 2 , and it assumes F A is a dipole 
with C = c 2 M 2 A = (1.1) (GeV/c) 2 . The reference model 
uses the derived Fp from Eq. [T3J and assumes that F v = 
F^ = at all Q 2 . We then define: 



IdQ 2 ^ 



JdQ 2 



JdQ 



, d<j r . 



(16) 



where a\ e ^ is the reference model for vgn — > £~p or its 
anti-neutrino analogue and ag is the model to be com- 
pared to the reference. 



A. Kinematic Limits 

The neutrino and anti-neutrino CCQE processes have 
kinematic limits in Q 2 which depend on the final state 
lepton mass, m. These limits are 
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(17) 



where s is the Mandelstam invariant representing total 
center of mass energy and E* t and p* t are the center of 
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FIG. 2. Our estimate in the lepton leg leading log approxima- 
tion of the fractional difference between the electron and muon 
neutrino total charged-current quasi-elastic cross-sections, A 
as defined in Eq. 1151 as a function of neutrino energy. The 
negative difference means that the electron neutrino cross- 
section is larger than the muon neutrino cross-section. 



mass final state lepton energy and momentum. E% can 
be expressed in terms of invariants as 



Et = 



s + m 



M 2 



2^ 



(18) 



Figure [T] shows the effect of the kinematic limits. Not 
surprisingly, the effect is very large near the threshold 
for the muon neutrino and anti-neutrino reaction. These 
effects are accounted for in the description of the quasi- 
elastic process in all commonly used neutrino generators. 
However, it is worth noting that the difference in Q 2 
spanned by the two reactions can lead to large effects in 
varying form factors that significantly affect either the 
small or large Q 2 parts of the cross-section. 



B. Radiative Corrections 

To calculate the effect of radiative corrections on the 
total quasi-elastic cross-section, we follow the approxi- 
mate approach of calculating the leading log correction 
to order log Q/m, where Q is the energy scale of the inter- 
action process [5(. This approach has a calculational ad- 
vantage in investigating the differences due to the lepton 
mass, m because the lepton leg leading log only involves 
sub-processes where photons attach to leptons. The key 
result from this approach is that the cross-section which 
allows for the presence of radiated photons, ulll is re- 
lated to the Born level cross-section, as, by 
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(19) 



where E\ is the center-of-mass frame lepton energy. 

In the case of elastic scattering, the relationship in as 
between Ei and the scattering angle, Of, simplifies the 
calculation because there is at most one z in the inte- 
grand for which the cross-section does not vanish for a 
particular lepton angle: 

z= [2E e (M + E v ) (to 2 + 2ME V ) - 2 cos 2 6 t EiE v 



x ym 4 + 4E 2 (M 2 - m 2 sin 2 6 e ) - 4m 2 M 2 - Ara 2 ME v 

I [to 4 + AE V (E v (to 2 cos 2 6 e + M 2 ) + m 2 M)) . (20) 

We then obtain the remaining cross-section by integrat- 
ing Eq. [19] over the final state lepton energy. Note that 
this procedure only gives a prescription for evaluating 
da(E Vt true) I dQtrue'i however, the radiation of real pho- 
tons means that the relationship between lepton energy 
and angle and E v and Q 2 in elastic scattering will no 
longer be valid. The effect of this distortion of the elastic 
kinematics will depend on the details of the experimental 
reconstruction and the neutrino flux seen by the experi- 
ment, so the effect must be evaluated in the context of a 
neutrino interaction generator and full simulation of the 
reconstruction for a given experiment. 

The difference of the effect on the total cross-sections 
as a function of neutrino energy is shown in Fig. [5] We 
estimate a difference of approximately 10% over the en- 
ergies of interest in oscillation experiments. The largest 
differences fractional differences in cross-sections are at 
high true Q 2 and low neutrino energies. The magni- 
tude of the lepton leg correction to the muon neutrino 
total cross-section is smaller, roughly 0.4 times this dif- 
ference, so the larger effect is on the electron neutrino 
cross-section. 

Our estimation of the effect is surprisingly large at the 
relevant energies for oscillation experiments. Some por- 
tion of this difference in the total cross-section in Fig. [2] 
may be canceled by diagrams missing from the leading 
log correction in the lepton leg, such as box diagrams in- 
volving 14^7 exchange between the leptonic legs and the 
initial or final state, which will also depend on the final 
state lepton mass [17| . We stress that this is only an ap- 
proximate treatment which should be confirmed in a full 
calculation implemented inside a generator, and to date 
radiative corrections are not included in the commonly 
used neutrino interaction generators^ 



C. Uncertainties in Fy, F v and F A 

As noted above, the vector form factors Fy and Fy are 
precisely measured in charged lepton scattering 18]; how- 
ever, the axial form factor is still uncertain because neu- 
trino experiments that measure it do not agree amongst 
themselves or with determinations in pion electroproduc- 
tion as discussed above. Therefore the axial form factor 
will dominate any differences in the electron and muon 
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FIG. 3. The change in the fractional difference of muon CCQE 
cross-section and electron CCQE when tua is changed from a 
reference value of 1.1 GeV/c 2 in a range generously consistent 
with current experimental data. 

cross-sections due to uncertainties in leading form fac- 
tors. 

Figure [3] illustrates the change in the fractional differ- 
ence of muon and electron neutrino CCQE cross-sections 
when the axial form factor is varied by changing the as- 
sumed dipole mass in a range consistent with experimen- 
tal measurements. The size of the effect is of order 1% at 
very low energy and drops with increasing energy. This 
difference in cross-section may be accounted for in vari- 
ations of the axial form factor within the analysis of an 
experiment using a modern neutrino interaction genera- 
tor. 



D. Pseudoscalar Form Factor 

At low Q 2 , the pseudoscalar form factor does have 
a significant contribution to the muon neutrino CCQE 
cross-section, of nearly the same order of the leading 
terms. However, Eq. 1131 shows that the contribution will 
be suppressed for Q 2 ~ M 2 , and all terms involving Fp 
are suppressed by m/M and so the contribution to the 
cross-section is negligible for electron neutrinos. At low 
neutrino energies, the pseudoscalar form factor effect on 
the cross-section difference, A(E l/ ) is nearly as large as 
that of the kinematic limits. The effect of the form factor 
as a function of neutrino energy and Q 2 is different for 
neutrinos and anti-neutrinos. 

Current neutrino interaction generators include 
the effect of Fp shown in Eq. [13] under the assumptions of 
PCAC and that the Goldberger-Treiman relation holds 
for all Q 2 . Experimental tests of the Goldberger-Treiman 
relation have identified small discrepancies which imply 
that the left hand side of Eq. [12] is between 1% 
6% less than the right-hand side 
models sugg ests that this effect is likely to disappear at 
high <5 2 [3fj]. We examine the effect of varying Fp(0) 



and 

Guidance from 



by 3% of itself as a reasonable approximation to the 
possible difference due to this effect. A more signifi- 
cant difference may arise due to violations of PCAC. 
This has been directly checked in pion electroproduc- 
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FIG. 4. The effect of variations of Fp from the reference 
model which assumes PCAC and the Goldberger-Treiman re- 
lation. The plots illustrate the change cross-section differ- 
ence, A(j5„), between a varied model and the reference model. 
Possible violations of the G-T relation produce a negligibly 
small effect, even at low energy. The range of violations from 
PCAC allowed by current data would allow significantly larger 
changes. The effect of setting Fp to zero is shown for com- 
parison. 

tion studies 24] which can directly measure Fp(Q 2 ) in the 
range of 0.05 to 0.2 GeV/c 2 . Uncertainties in this data 
limit the reasonable range of pole masses in Eq. [II] to be 
between 0.6M ff and 1.5M„. Effects due to these possible 
deviations from PCAC and the Goldberger-Treiman rela- 
tion are shown in Fig. [I] along with the effect of assuming 
Fp = for comparison. 



E. Second Class Currents 

As noted in the introductory material, non-zero sec- 
ond class currents violate a number of symmetries and 
hypotheses, and are therefore normally assumed to be 
zero in analysis of neutrino reaction data and in neutrino 
interaction generators. For this study, we take a data 
driven approach and look at the effect of the largest pos- 
sible second-class current form factors, F v and F\ that 
do not violate constraints from this data. 

Vector second-class currents enter the cross-sections 
for neutrino quasi-elastic scattering always suppressed 
by m/M and therefore only appear practically in muon 
neutrino scattering cross-sections. Both vector and axial 
vector form factors give large contributions to the B(Q 2 ) 
term given in Eqs. 3] and [7] and therefore typically have 
very different effects, often even different in sign, for neu- 
trino and anti-neutrino scattering. 

The vector second-class currents are difficult to de- 
tect in most weak processes involving electrons because 
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FIG. 5. S(E„,Q 2 ), defined in Eq.[H as a function of Q 2 for 
several selected E v - The difference between including and not 
including the maximum allowed second class vector current 
("SCCy"), F V (Q 2 ) = 4.4f^(Q 2 ), is shown. 
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FIG. 6. Changes in the difference between the muon and 
electron neutrino cross sections due to including F v . 
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Constraint 

FIG. 7. A survey of constraints on the ratio F a (0)/Fa(0) 
with their uncertainties from: (1) Wilkinson's data 
compilation 38] , (2) the same Wilkinson compilation with 
a correction for short-range effects [381], (3) the method of 
Wilkinson applied only to the A — 20 KDR parameter s 1371 . 
[3S|, (4) ibid, with a correction for short-range effects [37l. l38j. 
(5) a derived limit from A — 12 beta decays [3^] and (6) a 
derived limit from A = 20 beta decays [37T]. The value used 
for F% (0) /Fa (0) in this study is shown by the dashed line. 

the process is generally suppressed by powers of m e /M. 
Therefore even very precise beta decay measurements 
have difficulty limiting the size of F v (0) to less than 
several times the magnitude of the regular vector form 
factors '31]. The best limits from beta decays currently 
limit F v (0)/F v (0) to be (0.0011 ± 0.0013)^ 2.0 ± 



2.4 32[. Studies of muon capture on nuclei can pro- 
vide modestly better limits, but at the expense of as- 
suming there are no axial second class currents (31). An 
analysis of anti-muon neutrino quasi-elastic scattering 
has been used to place limits of similar strength, but 
again under the assumption of no axial second class cur- 
rents and with an assumed Q 2 dependence, Fy(Q 2 ) = 
F v (fi)/{Q 2 + M 2 v ) 2 with a fixed M 3V of 1.0 GeV/c 2 ^. 
From the preponderance of the data, we choose to pa- 
rameterize the maximum size of the allowed vector sec- 
ond class current as F V (Q 2 ) — 4.4Fy(Q 2 ), which is not 
excluded by the results of any of the above studies. The 
effect of this is significant, particularly at low neutrino en- 
ergies and is shown in Figs. [5] and [6l Recall that the effect 
on the electron neutrino cross-section from F v is negli- 
gible, so this effect occurs almost entirely in the muon 
neutrino cross-section. 

By contrast, the axial second class current at zero 
Q 2 is reasonably well constrained by studies of beta de- 
cay. We derive our limits in the framework of the KDR 
parameters [IJ where there is a wealth of experimental 
data to constrain these parameters (35T - l38j and therefore 
derive a limit on F\(0). Figure [7] shows these experimen- 
tal constraints and the effect we allow in this study. 

We assume a dipole form for the variation in Q 2 
of the axial second class current as well, so that 
for the maximum allowed variation F\{Q 2 ) / Fa{Q 2 ) — 
F|(0)/Fa(0) = 0.15. Figure U shows the effect of includ- 
ing this allowed axial second class current on both the 
difference of electron and muon neutrino cross-sections 
and on the muon neutrino cross-section itself. It is signif- 
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FIG. 8. Top: Changes in the difference between the muon and 
electron neutrino cross sections due to including F%; Bottom: 
the change in muon neutrino cross-sections due to including 



icantly smaller than the effect of the vector second class 
current because the limits on these currents are more 
stringent. 
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IV. CONCLUSIONS 

Large differences between the electron and muon neu- 
trino quasi-elastic cross-sections exist at low neutrino en- 
ergies from the presence of different kinematic limits due 
to the final state lepton mass and due to the presence 
of the pseudoscalar form factor, Fp, derived from PC AC 
and the Goldberger-Treiman relation. These differences 
are typically accounted for in modern neutrino interac- 
tion generators. 

There are also significant differences due to radiative 
corrections, particularly in diagrams that involve photon 
radiation attached to the outgoing lepton leg which are 
proportional to logQ/m. These differences are calcula- 
ble, but are typically not included in neutrino interaction 
generators employed by neutrino oscillation experiments. 
If our estimate of these differences, of order 10%, is con- 
firmed by more complete analyses, then this is a cor- 
rection that needs to be included as it is comparable to 
the size of current systematic uncertainties at accelerator 
experiments @, Q . 

Modifications of the assumed Fp from PCAC and the 
Goldberger-Treiman relation and the effect of the form 
factors Fy and F\ corresponding to second class vector 



FIG. 9. Top and Middle: For the form factors not well con- 
strained and not accounted for in neutrino generators, a sum- 
mary of the magnitude of the fractional size of differences in 
the total charged-current quasi-elastic cross-sections between 
electron and muon neutrinos and anti-neutrinos as a func- 
tion of neutrino energy. For Fp the average of the magnitude 
of the PCAC violating effects are summed linearly with the 
magnitude of the Goldberger-Treiman violation effect. Bot- 
tom: The magnitude of the difference between v and v of the 
fractional differences which illustrates the size of apparent CP 
violating asymmetries in oscillation experiments. 



and axial currents, respectively, are not included in neu- 
trino interaction generators. A summary of the possible 
size of these effects, as we have estimated them, is shown 
in Fig. M 

These differences, particularly from the second class 
vector currents, may be significant for current [2|-|J| and 
future [39[ neutrino oscillation experiments which seek 
precision measurements of — > v e and its anti-neutrino 
counterpart at low neutrino energies. Previous work[33j 
has demonstrated sensitivity to these second class cur- 
rents in neutrino and anti-neutrino quasi-elastic muon 
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neutrino scattering, and future work with more recent 
data[13, HH and newly analyzed data[4(| may help to 
further limit uncertainties on possible second class cur- 
rents. 
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